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The effect of blocking between different species occurring in one dimension is investigated here 
numerically in the case of particles following branching and annihilating random walk with two 
offsprings. It is shown that in two dimension simulations confirm the field theoretical results with 
logarithmic corrections. In one dimension however if particles exhibit hard core interaction I confirm 
the very recent predictions of Kwon et al. [PRL 85, 1682 (2000)] that there are two different 
universality classes depending on the spatial symmetry of the offspring production characterized 
by 13 s ~ 0.5 and Pa ~ 2. Elaborate analysis of simulation data shows that the order parameter 
exponent /3 does not depend on initial conditions or on diffusion rates of species but strong correction 
to scaling is observed. By systematic numerical simulations the critical point properties have been 
explored and initial condition dependence of the dynamical exponents Z and a is shown. In the 
case of a random initial state the particle-density decay at the critical point follows the t~^^'^ law 
with logarithmic corrections. 



I. INTRODUCTION 

Branching and annihilating random walk (BARW) 
process one of the simplest prototypes of a particle sys- 
tem exhibiting nonequilibrium phase transitions. For 
small creation/annihilation ratio the system evolves to 
a completely empty state from vi^hich no return is pos- 
sible. This so called "absorbing state" lacks any fluctu- 
ations. For higher creation/annihilation rates, a steady 
state with some finite concentration of particles emerges 
with a phase transition. This phase transition has been 
investigated first by numerical simulations and later 
with field theoretical methods Q. In general, we define 
annihilation of k particles kA — > with the rate A and 
generation of m offsprings A {m + 1)A with rate a. 

In D > 2 dimensions the phase transition of these sys- 
tems can be described by mean-field theory but in low di- 
mensions the effect of fluctuations introduces more com- 
plex behavior. The k = 1,2; m — 1 case is the so-called 
directed percolation process; the continuous phase transi- 
tion of it defines a robust universality class. According to 
the hypothesis of Janssen and Grassberger [^|J|] all con- 
tinuous phase transitions to a single absorbing state in 
homogeneous systems with short-ranged interactions be- 
long to this class provided there is no additional symme- 
try and quenched randomness present. As it was shown 
by Cardy and Taiiber Q this is the generic class for all 
critical models with k = l,2\m — odd and indeed a large 
number of model systems [Q have been found to belong 
to this class. 

If fc = TO = 2 the parity of the number of particles is 
conserved locally under the reactions and another univer- 
sality class appears: the so-called parity-conserving (PC) 
or BARW2 class. Particles following BARW2 dynamics 
may also appear as kinks between ordered domains in 
systems exhibiting two absorbing states p],p|-pl| . For re- 
cent reviews see For k > 2 the field theory does 



not predict any new universality classes. 

Our knowledge of multiparticle systems is quite incom- 
plete yet. A very recent paper of Janssen shows that 
this usually leads to asymmetries between pairs of species 
and unidirectionally coupled systems arise. These sys- 
tems have been shown to describe certain surface rough- 
ening processes |l^Jl^, for example. 

The bosonic field theoretical work of Q| predicts a 
new universality class for BARW2 process of symmet- 
ric N species (N-BARW2), exhibiting parity conservation 
(to = fc = 2) of particles. They found that offspring cre- 
ations of the same type like A AAA are not relevant 
processes because new particles are annihilated immedi- 
ately. In the case of branching processes like A ABB 
and A BBA, all > 1 component models have the 
same Feynman diagrams as that of the A^ — > oo model 
and the critical point for Z) < 2 happens to be at ct = 
zero branching rate with the following critical exponents: 

v\\=2/d, v^ = l/d, Z = 2, a^d/2, /3 = 1. (1) 

The critical exponents are defined by the usual way. The 
singular behavior of the spatial and temporal correlation 
lengths are characterized by 

^cxa-"- , (2) 

rocCT-"!! . (3) 

The anisotropy or dynamical exponent is Z = v\\lv^, 
while the particle density near the critical point obeys 
the scaling law 

p(t,a)cxt-Xto''ii) , (4) 

such that in the t — > oo steady state limit the order pa- 
rameter scales as p oc ct^; hence 
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j3 = avii = aZv^ 



(5) 



A. Zero branching rate case {a — 0) 



scaling relation holds. At the critical dimension [Dc = 
2) , logarithmic corrections to density decay are expected. 
For D > 2 the critical point is also at a = but governed 
by mean- field exponents [Q. 

In this work I investigated, numerically, the critical 
behavior of the 2-BARW2 model in one and two dimen- 
sions with the motivation that the existing field theories 
of reaction-diffusion models deal with a coarse grained, 
bosonic description of the particle system and therefore 
particle exclusion has not been taken into account. In 
a recent paper we have shown that in the case of 
one-dimensional (ID) annihilating random walk of two 
species, the blocking between different species can cause 
marginal perturbations. Therefore it is an open ques- 
tion whether this can also be observed if we add particle 
creation to this process. 

In a more recent paper Kwon et al. p^ ] argue based on 
local mean-field approximation that in this case new uni- 
versality classes emerge indeed. I confirm the predictions 
of for the order parameter exponent by precise nu- 
merical simulations and show that this exponent does not 
depend on initial-conditions or diffusion ratios of species. 
Furthermore I investigate the very nontrivial time depen- 
dence of this model and show numerical evidence for the 
initial condition dependence of the dynamical exponents 
and logarithmic corrections. 



II. DEFINITION OF THE MODEL 

A random sequential processes have been implemented 
in D = 1 and D ~ 2 dimensional lattices possessing pe- 
riodic boundary conditions with one time step described 
as follows: 

(1) A particle is selected randomly 

(2) A nearest neighbor site is selected randomly 

(3) With probability 1 — cr the particle is attempted to 
move to the new site if that is empty 

(4) If the new site is filled with a same kind of particle, 
both of them are removed 

(5) With probability a two off-springs of different types 
from the ancestor are placed at two free nearest 
neighbor sites selected randomly. 

In this way I investigate the same model as The 
field theoretical model of jj] also deals with the same 
particle reactions albeit without exclusion. 



III. SIMULATIONS FROM RANDOM INITIAL 
STATE IN TWO DIMENSIONS 

The initial conditions for systems with linear sizes L 
have been set up in such a way that L^/2 randomly se- 
lected A or B particles are placed at random sites. 



First the field theoretical expectations [p^P] for the 
density decay {p oc \n{t)/t) were tested. To see the log- 
arithmic corrections, I plotted p{t) x t versus ln(t). As 
Fig. |l| shows the curves corresponding to different sizes 
show linear behavior for short times. Later the effect of 
finite system sizes causes this behavior to break down. 
The crossover happens for larger and larger times as L 
goes to oo. 
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FIG. 1. Density decay (in 2D) of the total particle number 
for different system sizes: L = 16, 32, 64, 128, 512 (from top 
to bottom). Throughout the paper t is measured in units of 
Monte Carlo sweeps (MCS) of the lattices. 

For the largest system investigated here (L = 512) a 
logarithmic fitting was applied in the 10 < i < 10* re- 
gion, resulting in 

pit) = (0.442(2) + 0.296(1) ln(t))/t (6) 

Hence the field theoretical prediction a = 1 with logarith- 
mic corrections has been confirmed by these simulations. 
The amplitude of the logarithmic term however differs 
from that of the pure annihilation model: A2 = I/Stt 
p9| . One may speculate that in the two-component sys- 
tem this value is exactly I/tt. 

B. Off-critical case (a > 0) 

Traditional, direct simulations of off-critical exponents 
usually produce less precise estimates than cluster sim- 
ulations for example, but now I could eliminate the two 
most important drawbacks of this method. The critical 
point is known to be exactly at cr = 0; therefore this 
does not cause uncertainty in the estimation of the order 
parameter exponent. Furthermore I use precise analysis 
of the data by calculating local slopes that enables us to 
follow corrections to scaling. 

The simulations with branching were performed in sys- 
tems of linear sizes: L = 512 in the neighborhood of the 
critical point (ctc = 0) at ct = 0.0075, 0.01. ..0.17 (Fig. |). 
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About 100 samples were run from independent random 
initial states and averaged over for each a. 
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FIG. 2. Density decay (in 2_D) of the total particle number 
for a = 0.0075, 0.01. ..0.17 (from bottom to top) in a system 
of size L = 512. 

From the results at the critical point (Fig. |^) one can 
see that this system size is large enough for finite size 
corrections to be negligible for t < MCS investigated 
here. That means we can assume that in my steady state 
simulations the correlation length does not exceed L. The 
particle density curves were averaged following saturation 
and the effective pef / (c) exponents were calculated as 



Pefficr) = 



\npi - lnpi_i 
lufT," — lncr,_i 



(7) 



providing an estimate for the true asymptotic behavior 
of the order parameter 



cr^O 



(8) 



By applying a linear extrapolation as in Fig. |3| a very 
good agreement with the field theoretical prediction for 
the order parameter exponent f3 — 1.003(3) can be found. 
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FIG. 3. Effective order parameter exponent (/3) results in 
2D systems of linear sizes L = 512. 



IV. SIMULATIONS IN ONE DIMENSION 

A. The zero branching rate case (a — 0) 

In [0 we showed that in the case of special initial 
conditions 

A. .A. . .B.B. . .B. .B.A A B. .B. . . 

that assures that the time evolution never freezes (i.e. ev- 
ery particle has a neighboring particle of the same type 
to annihilate with) the density decays non-universally in 
an initial-condition dependent way. In this case parti- 
cles can be regarded as boundaries of compact domains 
of different species that cannot overlap each other. This 
happens in the genaralized Domany-Kinzel cellular au- 
tomaton, for example, [2^ . 

The decay depends on the initial density (and hence on 
the initial interparticle separation length) of the species. 
We argued on the basis of symmetry considerations and 
by using an analogy to the works of Kaiser and Turban 
pi| that this decay follows the powerlaw with continu- 
ously changing exponents. We also showed that in the 
case of asymmetric initial conditions, where a small seed 
of B particles is inserted in the see of A's, the cluster- 
survival probability of B's (characterized by the expo- 
nent 5) depends continuously on the initial density of 
A's similar to the case of the pair contact process [ p6| . 
Therefore in this work I investigated whether the static 
exponents show any initial-condition dependence as well 
and the form of the scaling law that connects them to 
the dynamic ones. 

In the case of random initial conditions particles segre- 
gate into coarsening mosaic of alternating domains con- 
taining alternating sequences of A's and B's. By map- 
ping this onto diffusion-controlled two-species annihila- 
tion {A + B ^ 9) Krapivsky and Redner argue that 
the particle density decays as i^^/'* |p3| . 

I tested this by simulating the process on L = 10^ lat- 
tices up to i = 4 X 10^ MCS. The local-slopes curves of 
the density decay defined as 



■ln[p{t)/pit/m)] 
ln(m) 



(9) 



(where usually I use to = 8) converge to ^ 0.25 indeed 
(Fig. ^), but again like in jl^ dependence on the ini- 
tial conditions can be observed. Now however, the local 
slopes do not saturate in the i — > oo limit as in the case 
of "pairwise" initial conditions ||l^ but show logarith- 
mically decreasing behavior as t — > oo. This is not an 
artifact of the finite system sizes since the same behavior 
can be seen on much smaller {L = 10^) lattices too. The 
origin of the logarithmic correction to scaling is not clear, 
but it is likely that marginal perturbation of species on 
each other causes it. Note that logarithmic correction to 
time dependence of the interparticle distances in case of 
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the A + B ^ model have also been shown by simula- 
tions and explained to be the consequence of density 
inhomogeneity inside the domains. 
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FIG. 4. Effective exponent of the particle decay in the one 
dimensional BARW2 with two species. The initial condition 
is randomly distributed A's and B's along the system of size 
L — 10^ sites. The solid line corresponds to initial densities: 
Pa ~ Pb ~ 0.375, the dashed line to: pA = Pb = 0.25, the 
dotted line to pA = pB ~ 0.125. The dot-dashed line corre- 
sponds to L = 10* and pA = Pb = 0.25, while the long-dashed 
line corresponds to L = 10* and pA ~ Pb = 0.125 

This anomalous decay behavior cannot be described 
by the conventional field theoretical description of a par- 
ticle system ^] in D = 1 that omits site restrictions and 
predicts a = 1/2. Also the paper of Kwon et al. p^ ] 
based on local mean-field approximation does not show 
any initial-condition dependency but predicts simply the 
a = 1/2 result. 

The dynamical exponent Z has been determined 
through the measurement of the characteristic time t{L) 
at the critical point. In this study I define r(L) as the 
time it takes for some portion of the sample to enter the 
absorbing state. In the case of pairwise initial conditions 
this portion was 50% but in case of random initial condi- 
tions, where the evolution is very slow in finite systems, 
I chose this portion to be 10% in order to get results for 
L > 1000 sizes. We expect the following finite size scaling 
at cr = 



t(L) oc 



(10) 



The data points have been determined within the L = 
16 — 4096 size region and to make precise extrapolation, 
local-slopes analysis has been performed. In Fig.|| I have 
plotted 



In-; 



In-; 



In Li — In L, 



(11) 



as a function of 1/L. As one can see, in the pairwise 
initial condition case Z^fj converges to Z — 1.82(2) as 
1/L — » quite smoothly. However in the random ini- 
tial condition case a huge correction to scaling can be 



observed and Ze/j seems to extrapolate to ~ 4.0(2), 
which is in accordance with the density decay law. 
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FIG. 5. Zeff vs 1/L at (T — 0. The square and 
star symbols correspond to pairwise initial conditions with 
pA ~ pB ~ 0.125, 0.25 respectively. The circles correspond to 
random initial conditions with pA ~ pB ~ 0.0625. 

The Z — 1.82(2) value is in agreement with our previ- 
ous finding where we showed that in the case of pair- 
wise initial conditions the density decays as ^-0.545(10)^ 
We argued there that the interparticle distance of AA 
and BB pairs Iaa = Ibb cx: p(t)~^ increases with the 
same power-law as the region of confinement that can be 
characterized with the exponent 1/Z. Now I could not 
obtain significant initial-condition dependence in Z ow- 
ing to the lack of such high precision simulations that we 
achieved in ||l^. Note that the usual scaling relation be- 
tween the single cluster spreading exponent z (describing 
the mean distance from the origin R{t) cx i^/2) and the 
dynamical exponent Z — 2/z is violated. This may be 
understood if we realize that in the seed growing process 
there is no symmetric, mutual exclusion between differ- 
ent species therefore the exponent z cannot describe the 
multisurface effect that arises in the case of pairwise ini- 
tial conditions. Neither of these Z values is in agreement 
with the assumptions of jl^ and the field theoretical re- 
sults for a ID system without particle exclusion M. 



B. The off-critical case (a > 0) with symmetric 
(static) branching 

The effect of the spatial symmetry in simple one- 
dimensional BARW processes was first investigated in 
psf and it was shown that reflection-symmetric (static) 
branching causes reentrant phase diagrams unlike the 
asymmetric (dynamic) branching process. First with the 
same special, pairwise initial conditions as in our recent 
paper , I performed steady state simulations such that 
offsprings are generated symmetrically on the left and 
the right sides of parents A — > BAB, B — > ABA (static 
branching). Kwon et al. [|l8| reported different critical 
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behavior than what was expected by applying the results 
of the bosonic field theory Q for this particle system. 

Similar to the two-dimensional case, the density decay 
was followed carefully and averaging was done within a 
long time period, following the steady state having been 
built up. Figure ^ shows my simulation results of the or- 
der parameter /3 in cases of different diffusion rates. As 
one can see, all the curves show strong (nonmonotonic) 
correction to scaling in the neighborhood of cr = but 
they all tend to the value (3 = 0.5. The overshooting 
of the effective exponent signals possible logarithmic cor- 
rection to scaling. The relaxation times to steady state 
were t < 10^ MCS in these simulations. 

I found that in the neighborhood of the critical point 
the steady state concentration is unchanged if I use larger 
lattice sizes (L = 4 x 10^ instead of 10^), therefore the 
strange shape of the curves cannot be an artifact of the 
finite system size of the simulations. Note, that had we 
done a simple fitting onto the data points we could have 
estimated (i ~ 0.55 instead of /3 = 0.5. 
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case, a'" ~ 0.25 ~ ~ 0.5/2 for random initial con- 

ditions. 

The scaling form (jj) can be checked by plotting pf^ 
against ta'^w . As Fig. ^ shows, one can obtain a fairly 
good data collapse with v^^^ = 2 and = 0.25 (corre- 
sponding to random initial conditions) for large times. 
For small times the collapse is less good owing to the 
effect of logarithmic corrections to . 




10"' 10"' , 10° 

to" 



FIG. 7. Data collapse for off-critical simulations of the sym- 
metric branching model according to the scaling form (^) with 
— 1/4 and i/^W = 2. The different curves correspond to 
a = 0.0005, 0.001, 0.002, 0.005 and random initial conditions. 



For pairwise initial conditions a very good data col- 
lapse can be achieved with the standard scaling form (^) 
using aP = 0.545, vF, = 0.915 



FIG. 6. Effective order parameter exponent (/3) results for 
the symmetric branching case. Circles correspond to the equal 
diffusion case, squares and triangles to dA = cis/2 nonequal 
diffusion rate case of species. 



C. The off-critical case {a > 0) with asymmetric 
(dynamical) branching 



The results are plotted as functions of ct°'^ determined 
experimentally such that the curves approach cr = in 
an asymptotically straight line manner. This means that 
the leading correction to the scaling exponent is ^ 0.5. 
As one can see, the curves converge to /? — 0.50(1) in 
agreement with the predictions of pq ]. The simulations 
with other initial densities and random initial distribu- 
tions of particles resulted in the same steady state values 
and hence f3 is unchanged. 

Comparing the Z's of the previous section with the 
= 1/2 of |l^] determined by static finite size scal- 
ing for the symmetric branching model we predict for 
pairwise initial conditions i/^ = Ziy± = 0.915(2) and 
i/^ = 2.0(1) for random initial conditions. 

The scaling law of density is fulfilled within error limits 
aP ~ 0.545 ~ P/'^f'i = 0.5/0.915 for the pairwise initial 



For the same special, pairwise initial conditions as in 
our recent paper [ pTf I performed steady state simula- 
tions such that offsprings were generated on the left or 
right sides of parents A ABB, A BBA, B BAA 
and B — > AAB. The relaxation time in this case was 
very long, typically t > 10^ MCS that prohibited getting 
closer to the critical point than a = 0.03. I performed 
these simulations in system sizes of L = 10^ sites and 
averaged over 100 samples. 

As one can see in Fig.|| the (3eff extrapolates to 
2.05(10) linearly in agreement with what was deduced 
from the local mean-field approximation by Kwon et al. 
p8[ . This value however, differs from the field theoretical 
predictions (/3 = 1) [Q obtained for the coarse grained, 
bosonic description of the particle system, which neglects 
hard-core interactions. 
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FIG. 8. Effective order parameter exponent results for 
tlie asymmetric brancliing case. 

I investigated the decay by starting from other different 
initial conditions as well. The steady state concentration 
and hence the exponent (3 was found to be the same, in- 
dependent of the initial density of particles. In the case 
of random initial distribution of particles the relaxation 
was found to be even longer reflecting the fact that ini- 
tial alternating sequences slow down the evolution to the 
steady state but the final concentration was the same. 
Again by comparing the Z's of Sec. IV A with v± = 2 oi 



mm determined by static finitesize scaling for the asym- 
metric branching process, we can predict for pairwise ini- 



tial conditions 



p 



Zv±_ = 3.66(2) andj/| = 8.0(4) for 
As in Sec. IV B the standard 



random initial conditions, 
scaling relation of the density is fulfilled for both initial 
cases within error limits: ~ 0.545 ~ ~ 2/3.66 

~ 0.25 ~ ~ ^/^ ^^"^ scaling law can be 

confirmed by data collapse (see Fig. ^). 
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FIG. 9. Data collapse for off-critical simulations of the sym- 
metric branching model according to the scaling form (^) with 
= 0.545 and t'j'J = 3.66. The different curves correspond 
to a = 0.03, 0.035, 0.04, 0.05 and pairwise initial conditions. 

In the case of dynamic branching Kwon et al. p!^ ] 
claim a slight dependence (a few percent) of the static 



exponents on the diffusion rates of different species. The 
relaxation time in this case is so huge that it hindered 
achieving such a precision by simulations. 



D. The off-critical case (a > 0) without exclusion 

Finally I tested the field theoretical results in one di- 
mension in such a way that I allowed particles of differ- 
ent types to exchange sites upon meeting. The steady 
state density has been determined for different values of 
cr around zero in lattices with L = 2 x 10"* size. As one 
can see in Fig.|l0| the middle curve (corresponding to this 
case) approaches to cr = linearly indeed in agreement 
with the (3=1 field theoretical result. This can be veri- 
fied by considering P^ff as well. 



V. CONCLUSIONS 

In this study I confirmed numerically the field theoret- 
ical predictions for the two dimensional branching and 
annihilating random walk process of two species. The 
logarithmic corrections to the power-law density decay 
function have been determined. 

In one dimension I gave numerical evidence that the 
bosonic renormalization group field theoretical predic- 
tions cannot be applied in contrast to widespread be- 
liefs for systems with exclusion and the order parameter 
exponent depends on the spatial symmetry of the off- 
spring production process. As Fig. |l^ shows, the steady 
state density is much lower if offsprings are created on 
the same side of the parent against the case when they 
are separated by the parent resulting in different scaling 
exponents j3 = 2 versus /? = 1/2. 
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FIG. 10. Steady state density in the one-dimenional 
2-BARW-2 model. Circles correspond to the asymmetric 
branching with /3 = 2, diamonds to the symmetric branch- 
ing with 13 = 1/2 and stars to the model without exclusion 
{(3 = 1). 



6 



This particle exclusion effect results in new universal- 
ity classes in one dimension as summarized in Table || and 
the static exponents are insensitive to the initial condi- 
tions. In the case of static branching the precise numer- 
ical analysis of the order parameter data does not show 
dependence on the diffusion ratios either. 

The density decay exponent starting from a random 
initial state (a) has been shown to agree with predictions 
of but initial density-dependent logarithmic correc- 
tions to scaling can be observed. The dynamical expo- 
nent Z has been found to depend on the initial conditions 
too. For pairwise initial conditions it is in agreement 
with our earlier study . All these results for dynami- 
cal scaling deviate from those obtained by the continuum 
field theory and assumed in [p^ . 

This study gives further evidence to the claims of Kwon 
et al. that in one dimension a series of new univer- 
sality classes emerge in nonequilibrium particle systems 
exhibiting hardcore interactions and continuous phase 
transitions. The relative simplicity of this model and 
the precision of the critical exponents obtained may pro- 
vide a good starting point to more elaborate analytical 
investigations of systems with hardcore exclusion. 

Acknowledgments 

The author would like to thank Z. Racz, S. Redner 
and U. Tauber for stimulating discussions. Support from 
Hungarian research fund OTKA (Nos. T-25286 and T- 
23552) and from Bolyai (No. BO/00142/99) is acknowl- 
edged. 



[1] P. Grassberger, F. Krause and T. von der Twer, J. Pliys. 

A:Math.Gen., 17, L105 (1984). 
[2] H. Takayasu and A. Yu. Tretyakov, Phys. Rev. Lett. 68, 

3060, (1992). 



[3] I. Jensen, Phys. Rev. E 50, 3623 (1994). 

[4] J. L. Cardy and U. C. Tauber, J. Stat. Phys. 90, 1 (1998). 

[5] H. K. Janssen, Z. Phys. B 42, 151 (1981). 

[6] P. Grassberger, Z. Phys. B 47, 365 (1982). 

[7] See examples in : J. Marro and R. Dickman, Nonequi- 
librium phase transitions in lattice models, Cambridge 
University Press, Cambridge, 1999. 

[8] N. Menyhard, J.Phys.A:Math.Gen., 27, 6139 (1994). 

[9] M. H. Kim and H. Park, Phys. Rev. Lett. 73, 2579, 
(1994). 

[10] H. Hinrichsen, Phys. Rev. E 55, 219 (1997). 

[11] K. E. Bassler and D. A. Browne, Phys. Rev. Lett. 77, 

4094 (1996). 
[12] H. Hinrichsen, Adv. Phys. 49, 815 (2000). 
[13] N. Menyhard and G. Odor, Brazilian J. of Physics 30, 

113 (2000). 

[14] H. K. Janssen, preprint, cond-mat/0006129. 

[15] U.C. Tauber, M.J. Howard, and H. Hinrichsen, Phys. 

Rev. Lett. 80, 2165 (1998); 
[16] H. Hinrichsen and G. Odor, Phys. Rev. Lett. 82, 1205 

(1999) . 

[17] G. Odor and N. Menyhard, Phys. Rev. E. 61, 6404 

(2000) . 

[18] S. Kwon, J. Lee and H. Park, Phys. Rev. Lett. 85, 1682 
(2000). 

[19] B.P. Lee, J. Phys. A 27, 2633 (1994). 

[20] H. Hinrichsen, Phys. Rev. E 55, 219 (1997). 

[21] C. Kaiser and L. Turban, J. Phys. A:Math.Gen 27, L579 

(1994); C. Kaiser and L. Turban, J. Phys. A:Math.Gen 

28, 351 (1995). 

[22] A. A. Ovchinnikov and Ya. B. Zel'dovich, Chem. Phys. 
28, 215 (1978); S. F. Burlatskii and A. A. Ovchinnikov, 
Russ. J. Phys. Chem. 52, 1635 (1978); D. Toussaint and 
F. Wilczek, J. Chem. Phys. 78, 2642 (1983); K. Kang and 
S. Redner, Phys. Rev. A 32, 435 (1985); M. Bramson and 
J.L. Lebowitz, Phys. Rev. Lett. 61, 2397 (1988). 

[23] P. L. Krapivsky and S. Redner, private communication. 

[24] F. Leyvraz and S. Redner, Phys. Rev. A 46, 3132 (1992); 

[25] S. Kwon and H. Park, Phys. Rev. E 52, 5955 (1995). 

[26] I. Jensen and R. Dickman, Phys. Rev. E 48, 1710 (1993); 
L Jensen, Phys. Rev. Lett. 70, 1465 (1993). 



process 


^1 




Z 


a 


P 


nonblocking 


2 


1 


2 


1/2 


1 


symmetric 


2.0(1)/0.915(2) 


0.50(1) 


4.0(2)/1.82(2)* 


0.25(1)/0.55(1)* 


0.50(1) 


asymmetric 


8.0(4)/3.66(2) 


2.0(1) 


4.0(2)/1.82(2)* 


0.25(1)/0.55(1)* 


2.05(10) 



TABLE L Summary of critical exponents in one dimension. The non-blocking data are quoted from [4]. The predictions for 
v± of the blocking models are cited from [18]. Data divided by "/" correspond to random vs. pairwise initial condition cases. 
Exponents denoted by * exhibit slight initial density dependence [17]. 
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